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Abstract 

The implications of the two-pocket Fermi surface for macroscopic quantum phenomena are considered. 
We demonstrate that in the case of the two-pocket Fermi surface the g-wave pairing is closely related to 
the d-wave one. As a result two macroscopic condensates arise. The Josephson tunneling for such two- 
component system has very special properties. We prove that the presence of the g-wave does not contradict 
the existing experimental data on tunneling. We also discuss the possible ways to experimentally reveal the 
g-wave component. 

1 Introduction 

There is a lot of controversy about the shape of the Fermi surface in cuprate superconductors. 
In the early days it was believed that it is a small Fermi surface of the doped Mott insulator 
[]]]. Later many of the results from photoelectron spectroscopy (PES) have been interpreted as 
favoring a large Fermi surface in agreement with Lattinger's theorem ||. On the other hand most 
recent PES data[[3], ^, ^, |6| once more give indications of a small Fermi surface for underdoped 
samples. In the present paper we consider the scenario with a small Fermi surface consisting 
of hole pockets around (±7r/2, ±7r/2), see Fig. la. It is widely believed that the t — J model 
describes the main details of the doped Mott insulator. To fit the experimental hole dispersion 
one needs to extend the model introducing additional hopping matrix elements t', t" (see. e.g. 
Refs.[[7|, H), but basically it is the t — J model. Superconducting pairing in the t — J model 
induced by spin- wave exchange has been considered in the papers [JTU, [TI| . It was demonstrated 



|TT| that there is an infinite set of solutions for the superconducting gap. All the solutions have 
nodes along the lines (1,±1), see Fig. la. (It is very convenient to use the magnetic Brillouin 
zone, but it can certainly be mapped to the full zone.) Using translation by the vector of inverse 
magnetic lattice the picture can be reduced to two hole pockets centered around the points 
(±7r/2, 7r/2), see Fig. la. The superconducting pairing is the strongest between particles from 
the same pocket, and the lowest energy solution for the superconducting gap has only one node 
line in each pocket. Having this solution in a single pocket one can generate two solutions in the 
whole Brillouin zone taking symmetric or antisymmetric combinations between pockets. The 
symmetric combination corresponds to the d-wave (Fig. lb), and the antisymmetric combination 
corresponds to the g-wave (Fig.lc) pairing. We would like to note that the possibility to generate 
new solutions taking different combinations between pockets was first demonstrated by Scalettar, 
Singh, and Zhang in the paper |I2 |. The energy splitting between the d- and g-wave solutions has 



been investigated numerically in the Ref.[^T|. The g-wave solution disappears and only d-wave 



one survives as soon as the hole dispersion is degenerate along the face of the magnetic Brillouin 
zone. Actually in this situation there are no pockets and one has a large open Fermi surface 
at an arbitrary small hole concentration. However for small well separated pockets the d- and 
g-wave solutions are almost degenerate. In pure t — J or t — t' — t" — J model the d-wave solution 
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always has the lower energy. However if one extends the model including nearest sites hole-hole 
Coulomb repulsion the situation can be inverted. The nearest sites repulsion does not influence 
the g-wave pairing and substantially suppresses the d-wave pairing. So it is quite possible that 
the real ground state has the g-wave superconducting gap. 



2 Ginzburg-Landau free energy 

In the present paper we consider the scenario with hole pockets well separated in k-space. Ac- 
cording to the described above microscopic picture we should consider simultaneously the d- and 
the g-wave pairing. In general terms such situation was considered a long time ago for conven- 



tional superconductors by Leggett fL3| . Let us formulate an effective Ginzburg-Landau (GL) 
theory. In the first approximation we can neglect the interaction between pockets and intro- 
duce two macroscopic condensates corresponding to the pockets, = \^>i\ e 4 ^ 1 , \I/ 2 — I ^2 1 e"^ 2 , 
where = |\l/ 2 | — |^| = is the number density of the condensate holes. The 

Ginzburg-Landau free energy of a such system in an external magnetic field can be written as 

F = /^((V-^A^I^— |(V-^A)^ 2 | 2 - (1) 
J V2m* IV he ' 1 2m* IV he ' 21 y ' 

a|^| 2 + o|^| 4 - a|* 2 | 2 + o|^ 2 | 4 + f - ] dV. 
A small interaction between pockets can be described by adding to the free energy ([!]) the term 

Fint = 1 J (*J*2 + dV 2 7^ |*i| |* 2 | cos(0! - 2 ), (2) 

where 7 is a small parameter of the interaction, 7 < a. The total bulk energy of the cuprate 
superconductor equals F v 9 = F + Fj n t, and the equilibrium values of the order parameters are 

|^| 2 = |# 2 | 2 = |^| 2 = ( a + |7|)/2& « a/26. (3) 

The ground state phase difference A0 = 2 — 0i is determined by the sign of 7: if 7 > 0, 
then A0 = tt ( g-wave ); if 7 < 0, then A0 = ( d-wave ). It is also convenient to introduce 
the d- and g-wave condensates ^ d = ^1 + ^ 2 = 2|^|cos(A0/2)e i<?i and ^ g = ^1 - ^ 2 = 
2|\E'|sin(A0/2)e l9i_l7r / 2 , where 20 = 0i+0 2 . So the ground state has either d- or g-wave symmetry: 
if 7 > then ^ d = 0, V g ^ 0, if 7 < then ^ g = 0, ^ d ^ 0. 



3 Tunnel junction 



Let us consider a Josephson tunnel contact of a conventional s-wave superconductor (the order 
parameter \l/ s = |^ s | exp(i<p s )) and a high-T c cuprate superconductor. Constant supercurrent 
j is maintained through the contact. The bulk free energy of the conventional superconductor 
equals 

B T 



h2 ^ 2ze . , T - 
(V--A)* 



2m 



a, m t 



and the surface free energy related to the contact is 



-Ai(n) + -A 2 (n) 



+ )dV, 

671 / 



dS, 



(4) 



(5) 



where n is a unit vector orthogonal to the surface of the contact and directed from conventional 
superconductor to cuprate, Ai(n) and A 2 (n) are the tunneling matrix elements from the first 
and the second pockets correspondingly. £ is an axis parallel to n, and ^2+, and \I/ S _ 
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are values of the condensates at the contact (£ = ±0). It is convenient to use following simple 
parameterization for the tunneling matrix elements 



Ai(n) 
A 2 (n) 



C{nl 
C(nl 



fiJ)e- a(Tta_n » )2 , 
n 2 )e~ a( ' nx+ny ^ 2 , 



(6) 



where C and a are some parameters depending on the tunneling probability, the shape of the 
Fermi surface, etc. The axes x and y are directed along crystal axes a and b of the cuprate. 
We stress that eq.(^) is just a parameterization having correct symmetry properties. Let us 
introduce also the d- and g-wave tunneling matrix elements: A<i(n) = Ai(n) + A2(n) and A ff (n) = 
Ai(n) — A2(n). The symmetries of these matrix elements shown at Fig.2 are the same as the 
symmetries of d- and g-wave gaps in the momentum space, Fig. lb,c. After variation of the 
total free energy F = F v 9 + F v + Fg with respect to *2> *s and vector potential A we find: 
the GL equations for conventional and high-T c superconductors 
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*i - a$>i + 2b *i + 7*2 = 0, 



(7) 



* 2 - a* 2 + 26 |*2 1 ^2 + 7*i = 0, 



the boundary conditions for these equations 
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(8) 



and the currents in the cuprate and conventional superconductors 
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Substituting the boundary conditions (|§|) into eqs. 
through the contact 



we find the expression for supercurrent 



j = ^(n)^*^ - *:_* 1+ ] + ^A a (n)[¥._n, 



*:_^ 



(10) 



We will consider the distances larger than the cuprate superconducting correlation length. At 
these distances the magnitudes |*i| and |*2| are equal to the constant given by eq.(^), and only 
the phases of condensates are £ dependent. The Josephson current through the contact can be 
written as 

4e 

j = |^| |* s | [A^n) sin(0 1+ - S _) + A 2 (n) sin(0 2 + - &-)], (H) 

where i+ and 2 + are values of the phases at the contact (£ = ±0). We consider the 
situation without an external magnetic field in the contact. The magnetic field caused by the 
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supercurrent is very small, therefore we set A = and reduce eqs.(|7]) for ^ and ^2 to ones for 
the phases 



h 2 d 2 A(j) . h 2 d(f) dA(j) 



+ ^ 



2m* d£ 2 m* d£ d£ 
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(12) 
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1 + j^j cos A(f)) = 0, 



with 



= (0i + 02)/2 and A0 = 4>2 — 4>i- The equilibrium value of A0 is A0o with A0o = 
at 7 < ( d-wave ), and A<fto = tc at 7 > ( g-wave ). If the current is small we can linearize 
eqs.(O) with respect to A0 — A(f> . This gives that 4> = (<pi + 4>2)/2 = const, and 



A0 



>i = A0 O + Aexp(-£/l 
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where / 7 = h/JA \^ \ m*. It is convenient to rewrite the supercurrent ( pTT| ) in the form 



4e 
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(13) 
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where 9 = (p — </> s _ and A0 C is the value of Acj) at the contact: A<p c = A for 7 < ( d-wave ) 
and A(p c = it + A for 7 > ( g-wave ). The solution we found depends on the two arbitrary 
constants 9 and A. If the supercurrent through the contact is fixed then the eq. (|TJ) gives one 
relation between these constants. Second relation follows from the minimum of total free energy. 
After substitution of the solution ( |i~3"D into expressions ([[]), (H), @ we find the part of the total 
free energy depending on 9 and A 



m 2 k^a 2 -2m 



\h cos 9 cos 



+ A sin 9 sin 



. A0 C 



(15) 



Here K~ 



Z 7 | 7 | = n^J\ 1 \/m*/2. 

Let us consider at first the contact plane to be parallel to the axis a or axis b of the cuprate 
(n = (±1,0) or n = (0, ±1)). In this case X g = and Ad 7^ 0. If the ground state pairing has 
d-wave symmetry (7 < 0) we find from eqs. (Tj) and ([15]) the following condition of minimum 
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A = 0. 



j=const \ 

Therefore A = 0, and 9 is determined by the usual relation 



j = j cd sin 0, j cd 



h 



(16) 



(17) 



If the ground state pairing has g-wave symmetry (7 > 0) the condition of minimum is 



dFn 



dA 



This gives 

and the supercurrent is 



j=const 
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(19) 



(20) 
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One has to remember that the above equations are derived assuming that A <C 1 and therefore 
they are not valid when 9 is approaching ix/2. So the picture is that deep inside the cuprate we 
have g-wave pairing, but in the layer of width Z 7 near the contact there is admixture of the d-wave 
component. If current is increasing the admixture is increasing too. So the current pumps the 
g-wave into the d-wave in the layer of width Z 7 . At critical current there is only the d-wave at 
the contact. Therefore the critical current in this case is exactly the same as one for the d-wave 
ground state: j c = j cd . 

At an arbitrary orientation of contact plane with respect to the crystal axes both A^ and X g 
do not vanish. Nevertheless the d-wave tunneling probability is always larger than the g-wave 
one: |Ad(n)| > |A s (n)|. It is obvious from eqs.(||). Let us rewrite the current ( fT4[) as 



2y/2e 



2 + \t + (\ 2 d -\*) 



cos 



■ sin 



— a 



(21) 



where sin a = v^2A s sin 4^ / \J\ d + \ 2 g + (A^ — A~|) cos A0 c , and let us introduce the notations 



Jcd 



Jcg 



keeping in mind that j cd > j cg - Let us assume first that the 
ground state of the cuprate has the d-wave gap (7 < 0). Then at any supercurrent the system 
is in the d-wave everywhere including the contact (A<p c = 0), and from eq.(|2~l"D we conclude 
that the critical current is j c = j cd . Consider now the situation when the ground state of the 
cuprate has the g-wave symmetry (7 > 0). Then at j < j cg the system remains in g-wave state 
everywhere including the contact, A0 C = ir. Increasing the current above j cg we start to pump 
the g-wave into the d-wave in the surface layer of the width l 1 . When the current reaches j cd , 
according to eq.(|2T|) the phase takes the value A0 C = 0, and this is the critical point. So at the 
critical point in the contact the g-wave is completely pumped into the d-wave, and the critical 
current is exactly the same as in the case of the d-wave pairing: j c = j cd . 



4 SQUID 



The g-wave pairing can be revealed in the SQUID. Consider at first the geometry of the SQUID 
experiment , results of which are interpreted as a very strong evidence in favor of the d-wave 
pairing. The setup is shown schematically at Fig.3. Faces of the superconducting corner are 
parallel to crystal axes a and b of the cuprate superconductor so that A 9 (n) = and Arf(n) have 
opposite signs at the tunnel contacts: A^i = — \&i = Xd- The indexes 1,2 numerate the contacts. 
Using eqs.(|l4]),(|l"5"l) the total current and the free energy can be written as 



4e 

3= ji+32= — |#||* S |A 
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Here 6^, A i: and A0 ci = A0 O + A { (i=l,2) are the quantities introduced in previous section. 
Quantization condition is standard 



e l = e + 
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27rm, 



(24) 



where $0 — nch/e is the quantum of magnetic flux, and m is an integer number. 

Let us consider at first the case 7 < which corresponds to the d-wave pairing in the cuprate 
ground state. Minimum of the free energy (^) at fixed supercurrent (^2|) is defined by the 
equations 



2K~\W + 



j=const 



A d |^ a |sin(7r$/$ ) 
2cosfl 



Ai = 0. 



(25) 
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The solution of these equations is A\ = A 2 = 0, so the system remains in the d-wave state 
everywhere including the contacts, and the current equals 



j = j cd sin 9, j cd = — \^^ s X d \ 

n 



sm 7T 



$ 



(26) 



The plot of j c d is represented at Fig. 4a and is in agreement with experimental data [13]. 

Now consider the case 7 > which corresponds to the g-wave pairing in the cuprate ground 

state. Minimum of the free energy ( p3f) at fixed supercurrent fl22|) is defined by the equations 



j=const 
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j=const 

There are two solutions of these equations: 

symmetric : A\ - 
antisymmetric : A% - 



A\ cos 9\ — A 2 cos 9 2 
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(27) 
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The supercurrents and free energies corresponding to symmetric and antisymmetric solutions 
are 
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/K 1 . The indexes (+) 



(-) denote symmetric and antisymmetric solutions 
correspondingly. The physical state corresponds to minimum of the free energy and therefore 
at — 1/4 + m < $/$o < 1/4 + m the SQUID is in the antisymmetric state (Ax = — A 2 ), and at 
1/4 + m < $/$o < 3/4 + m it is in the symmetric state (Ai = A 2 ), m is an integer. One has 
to remember that eqs. ( plf) and (|29|) are derived at Ai <C 1, and they are not valid when 9 is 
approaching 7r/2 for the symmetric solution and when 9 is approaching for the antisymmetric 
solution. So at 7 > inside the superconductor we have g-wave pairing, but in layers of width 
Z 7 near the contacts there is admixture of the d-wave component. If the current is increasing the 
admixture is increasing too. So the current pumps the g-wave into the d-wave state in the layer 
of width lj. At the critical current there is only the d-wave at the contact. The critical currents 
corresponding to symmetric and antisymmetric solutions can be easily found from eq. fl22|) if we 
substitute Ai = A 2 = n or A\ = —A 2 = it. 



J. 
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(-) 



8e 
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W s \ d sin vr 



W„\ d cos U 



$0 



(30) 



The real critical current is j cg = maxjj^, )} The plot of j cg as a function of magnetic flux 
is given at Fig.4b. It is interesting that it has a period $o/2 which could be naively interpreted 
as flux quantization with an effective charge 4e. 



G 



So we see that the dependence of the SQUID critical current on the magnetic flux is different 
for different ground states of the cuprate superconductor. We considered above the geometry 
with contact planes parallel to crystal axes a and b of the cuprate. If the bulk ground state of 
the cuprate has d-wave pairing (7 < 0) the system remains in the d-wave everywhere including 
the contacts and the critical current is proportional to |sin (7r$/$o) I- If the bulk ground state 
of the cuprate has g-wave pairing (7 > 0) then the current pumps the g-wave into the d-wave in 
the layers of the width Z 7 near the contacts. The critical current in this case has a very unusual 
dependence on the magnetic flux presented at Fig.4b. 

Let us consider now a more general geometry of SQUID contacts when X g 7^ 0. We assume that 
it is a straight angle rotated with respect to crystal axes of cuprate, therefore A^i = —Xd2 = Xd, 
A 9 i = A 9 2 = X g , see Fig. 5. We assume also that the angle of rotation is not very small, so that 
Ad ~ A r Consider at first the case of the d-wave ground state of the cuprate superconductor 
(7 < 0, A0o = 0). The free energy is given by the sum of eqs.(|15"|) corresponding to the two 
contacts. The same is valid for supercurrent (|T4|). Minimization of the free energy at fixed 
supercurrent with respect to variation of phases Ai gives in this case a unique solution. At 
Ai <C 1 the solution can be found analytically 



A 9 1 




COS (7T$/$ ) 


2K^ 


1^1 


cos 9 



A 1 = A 2 = v ' uy . (31) 



So at X g ytz external current pumps the d-wave into the g-wave in the layer of width Z 7 near 
the contacts. Keeping in mind that A\ = A2 = A we can write the supercurrent at arbitrary A 
in the form 

j = j c . cos(a + A/2) ■ sin 6, (32) 
with sin a = X g cos(7r J^) / ^JX 2 d sin 2 (7r J^-) + A 2 cos 2 (7r^-) , and the critical current given by 



J C = ll^l^sin^vr^ + A^cos^vr^). (33) 

The plot of the critical current is given at Fig. 4c. The admixture of the g-wave at the contact at 
critical current is defined by a. It is equal to 100% at $ = 0, and it is equal to at $ = $o/2. 

Now we consider the last situation: geometry with A^i = — Xd2 = Xd, X g i = X 9 2 = X g 7^ 
(see Fig. 5) and the g-wave ground state of the cuprate superconductor (7 > 0, A0 O = n). 
Minimization of the free energy shows that only one solution exists, which at small Ai has the 
form 



A l = A 2 = — ,' v ; u/ . (34) 

2K, t cos# V ; 

It is similar to the symmetric solution (|28| ) for the geometry with X g = 0. It is interesting that 



in this case there is no an analog of the antisymmetric solution (|28|) . So, as usual the external 

current pumps the g-wave into the d-wave at the contact surfaces. The supercurrent at an 
arbitrary Ai, but with A\ = A 2 = A can be written as 

j = j c ■ sin(a + A/2)- sin 9, (35) 



with exactly the same j c and a as those for the d-wave ground state, given by eqs. (|32|), (|33|). 
The admixture of the d-wave at the contact at critical current is defined by a. It is equal to 
100% at $ = $o/2 and it is equal at $ = 0. So the critical current through the SQUID is 
given by eq.(^) and it is independent of whether we have the d-wave or the g-wave pairing in 
ground state. The plot of the critical current is given at Fig.4c. 

The following question arises. In the case of contact planes parallel to the crystal axes a 
and b the dependences of critical current on flux are different for different ground states: the 
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dependence is given by Fig. 4a for the d-wave ground state, and it is given by Fig. 4b for the g- 
wave ground state. On the other hand as soon as the contact is rotated (Fig. 5) the dependence 
is the same for the d- and g-wave ground states, and it is given at Fig. 4c. What is the angle 
of rotation necessary to change the regime? The origin of the difference in regimes is in the fact 
that at X g = there are both symmetric and antisymmetric solutions, see eq. (P3|), but at finite 
X g the anisymmetric solution disappears, see eq. (|34|). One can prove that with increasing X g 
from zero the free energy corresponding to the antisymmetric solution is increasing very fast, 
and at X g /Xd > j3 c , 

*~2^W <<1 ' <36) 
it becomes very high. So the regime is changed when the rotation angle (3 is larger than f3 c . The 
value of f3 c is very small because it is proportional to the tunneling amplitude. 



5 Conclusions 

We have considered the scenario with the small Fermi surface consisting of the hole pockets. 
The picture can be relevant to underdoped cuprate superconductors. The small Fermi surface 
together with mechanism of the magnetic pairing results in the possibility of having both the 
d- and the g-wave pairing. Energy splitting between these states is small. The ground state 
symmetry is defined by the interplay between the magnetic pairing and Coulomb repulsion. We 
have formulated an effective Ginsburg-Landau theory describing this situation. 

We demonstarate that in the Josephson junction or the SQUID consisting of cuprate super- 
conductor and conventional superconductor the supercurrent pumps the d-wave into g-wave (or 
the g-wave into d-wave) in the thin layer near the contact. This pumping influences the SQUID 
interference picture. If the bulk ground state has d-wave symmetry and contact planes are par- 
allel to crystal axes a and b of the cuprate the dependence of the SQUID critical current on the 
magnetic flux is shown at Fig. 4a. This is a well known picture for the d-wave superconductor. 
If the bulk ground state has g-wave symmetry and contact planes are parallel to crystal axes 
a and b of the cuprate the dependence of the SQUID critical current on the magnetic flux is 
shown at Fig. 4b. It is very unusual and could be naively interpreted as flux quantization with 
an effective charge 4e. If the contact planes are rotated with respect to crystal axes a and b at 
the angle j3 > /3 C , then the SQUID interference picture is the same for the d- and g-wave pairing 
in the bulk ground state, and this picture is presented at Fig. 4c. The angle (3 C given by eq. ([36]) 
is very small, since it is proportional to the tunneling amplitude. 
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FIIGURE CAPTIONS 



Fig. 1. a. Fermi surface in magnetic Brillouin zone which is equivalent to the two-poket Fermi 
surface (dashed line), b. Symmetry of the d-wave pairing in momentum space, c. Symmetry of 
the g-wave pairing in momentum space. 

Fig. 2. a. The d-wave tunneling amplitude, b. The g-wave tunneling amplitude, n x and n y 
are crystal axes of cuprate. 

Fig. 3. Geometry of the SQUID experiment with contact planes parallel to the crystal axes. 

Fig. 4. Dependence of the SQUID critical current on magnetic flux. a. The d-wave bulk 
ground state and contact planes exactly parallel to the axes a and b of the cuprate. b. The 
g-wave bulk ground state and contact planes exactly parallel to the axes a and b of the cuprate. 
c. The d- or g-wave bulk ground state and superconducting corner rotated by a angle (3 > f3 c 
with respect to the crystal axes a and b. 

Fig. 5. Geometry of the SQUID experiment with the superconducting corner rotated by an 
angle j3 with respect to crystal axes. 
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